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The recent direct photon data, especially the large elliptic flow from AuAu collisions at y^sjvjv = 
200 GeV and PbPb collisions at ^/sNN =2.76 TeV has been investigated with hadronic data con- 
strained 3+ ID ideal hydrodynamical models. It turns out that those data support an earlier 
starting of the collective motion (ro < 0.6 fm/c) and a delayed formation of three-flavour QGP 
( TQGP ~ 2 fm/c) in heavy ion collisions. The observed large elliptic flow of direct photon is still 
underestimated. A more realistic model with viscosity, an earlier to and event-by-event fluctuation 
is expected to improve the work. 

PACS numbers: 



I. INTRODUCTION 

The recent direct photon data [l| from Pb+Pb col- 
lisions at ^Jsnn = 2.76 TeV shows a big enhancement 
at low transverse momentum (pt) region, similar to the 
PHENIX observationQ . Also recently, a large elliptic 
flow (U2) of direct photons, as large as that of hadrons, 
has been observed, both by PHENIXQ and preliminarily 
by ALICE The large elliptic flow of direct photons 
challenges our understanding of heavy ion collisions. 

To explain the large elliptic flow, some people con- 
sider additional sources of photon emission. For exam- 
ple, conformed anomaly as a source of soft photons [H 
can produce a sizable elliptic flow. However, this addi- 
tion source does not vanish at high pt region and is ruled 
out by the PHENIX observation that Raa = 1 at high 
Pt ■ The contribution of pre-equilibrium Glasma stage 
to photon emission is estimated in where a number of 
approximation must be made due to the current incom- 
plete knowledge of this stage, which claimed this source 
may generate elliptic flow but hard to calculate. In Q, 
initial fluctuation has been considered and many useful 
qualitative results have been obtained. 

Now we would like to check how well photon data and 
hadronic data can be simultaneously explained. There- 
fore the 3+ ID ideal hydrodynamical model Q is chosen 
for AuAu collisions at i/sivjv = 200 GeV and [l^l for 
PbPb collisions at y/JWN =2.76 TeV. In both cases, hy- 
dro parameters are well fixed to reproduce excellently the 
bulk hadrons such as rapidity distribution, pt spectra and 
elliptic flow. More extensive test on charm production 
and jet quenching, HBT and so on has been performed 
with those hydro models. The difference between the 
two models, i.e., initial time tq = 0.6 fm/c for AuAu 
and To = 0.35 fm/c for PbPb will tell us more than a 
single model chosen. Since data on bulk hadrons have 
been successfully explained, in this work we will focus on 
how well photon data can be explained. 



'Electronic address: |liufm@iopp. ccnu.edu.cn] 



Additionally, we would beneflt from photon data and 
ask a simple question, "When is QGP formed in heavy ion 
collisions?" Here the QGP refers to a 3- flavour (locally) 
thermalized partons system to meet the observation of 
the saturated strangeness suppression factor 7^ and free 
baryon number at midrapidity for central collisions jTlj. 

Conventionally, this QGP matter is assumed to be 
formed at the beginning of collective motion, i.e., the 
electromagnetic radiation from QGP are often calculated 
from the hydro initial time tq |12h14| . However, in heavy 
ion collisions at very low energies, no QGP formation but 
hydrodynamical model is still widely employed. In Lan- 
dau hydrodynamical model [isj , the hydro initial time is 
as early as just after the little bang, certainly QGP for- 
mation should be later since it takes time for the system 
to realize local equilibrium, both thermally and chemi- 
cally. Hydro model is so widely employed, because the 
conervation law and the general relation e = e{p) to close 
hydro problem, hold for many kinds of matter. But we 
should distinguish QGP formation time tqgp from hydro 
initial time tq. 

How to estimate QGP formation time tqgp? The 
well-known thermal photon emission rates count photons 
emitted due to the secondary interactions in a QGP or 
a hadronic gas [TgI . [TtI . So in this work, 3-flavour QGP 
emission rate |l7| will be employed, and the beginning 
time of counting varies instead of the fixed value hydro 
initial time tq. The appropriate beginning time for this 
counting to mimic the recent data will provide us an es- 
timation of QGP formation tqgp . 

The paper is organized as following: after a brief in- 
troduction of calculate approach, we present the results 
and the interpretation, then conclusion and discussion. 



II. APPROACH 

In our previous work[l3j. we successfully explained the 
Pt spectra of direct photons from AuAu collisions mea- 
sured by PHENIX with the hydrodynamical model cho- 
sen in this work, and made a detailed investigation on 
the competition among the main direct photon sources. 
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The collisions between hard partons and the bulk mat- 
ter, make a photon source called as jet-photon conver- 
sion(JPC) whose contribution increases with system size. 
At the mean time, those collisions reduce the contribu- 
tion of another source, fragmentation photons, or the 
higher order contribution of prompt photons. The in- 
crease and reduction cancel each other, so that Raa = 1 
and 1)2 = at large pt and coincide with the measure- 
ment. Thus, we can simplify the direct photon calcu- 
lation, considering neither JPC nor energy loss. There- 
fore, we consider only two main sources of direct photon 
production, prompt photons and thermal photons. The 
pre-equilibrium contribution to photon emission is not 
considered separately in the calculation. We will discuss 
the effect due to this approximation. 

The prompt photon production is calculated as next- 
to-leading order cold nuclear collisions: 



dy(Ppt 



= TAB{h)Y, J dXadXbGa{Xa,M^)Gb{xi„M^) 
ab 



X -S{s + i+ u)[^{ab^ ^ + X) 
TT dt 



(1) 
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where TAsib) is thickness function at given centrality, 
G{x,M'^) is nuclear parton distribution function with 
nuclear shadowing effect and EMC effect, and the rest 
is similar to the calculation in pp collisions, i.e. [l8| . 

The thermal photon calculation is based on the above- 
mentioned hydrodynamical description of each collision 
system. The solution of hydrodynamical equation will 
provide us the temperature T and flow velocity at 
each space-time point of the system. Thus, the pt spec- 
trum of thermal photons reads 



dN^ 
dyd'^pt 



d^xT{E*,T) 



(2) 



where T{E*,T) is the photon emission rate at tempera- 
ture T, which covers emission from 3-flavour QGP and 
later formed hadronic gas. E* = p^u^, p^ is the four- 
momentum of a photon in the lab frame. Conventionally, 
d'^x makes the space-time integral over the whole space 
and the whole emission history of the system, so the time 
integral begins at hydro initial time tq. In this paper, we 
will try different choices. The best choice may provide 
an estimation to the QGP formation time tqgp- 

The elliptic flow of thermal photons, quantified by the 
second harmonic coefficient 



/ d(t)Cos{2(j))d^N^/dy d^pt 
Jd4^d^W^Jdy¥pt 



(3) 



comes from the Lorentz boost from the local rest frame to 
the lab frame. In the former frame, photons are emitted 
isotropically. The connection of the two frames is flow 
velocity u^. So both its strength and asymmetry in the 
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Figure 1: (Color Online) Photons from PbPb collisions at 
y/SNN = 2.76 TeV, centrality 0-40%. Calculated results with 
four beginning times are presented as four kinds of curves. 
Full squares are ALICE data [l|, i|. For more details, see 
text, (a): Transverse momentum spectrum of direct photons, 
(b): Transverse momentum spectrum of thermal photons and 
prompt photons, (c) The elliptic flow of direct photons, (d) 
The elliptic flow of thermal photons. 



non-central collisions, are crucial for the elliptic flow of 
thermal photons (l9| . 

To compare with the elliptic flow data, we should take 
into account of prompt photons, though they have V2 
zero. Thus the elliptic flow of direct photons reads 



TvJipt), 



where [vF ) refers to dN^ / dyd^pt {dN^ / dyd^pt) 



(4) 



III. RESULTS 

In Fig. 1 is shown the midrapidity photons from PbPb 
collisions at ./sWn = 2.76 TeV of 0-40% centrality Full 
squares are ALICE data Four curves present re- 

sults from four beginning times, r =0.35 fm/c (solid line), 
0.75 fm/c (dashed line), 1.15 fm/c (dotted dashed line) 
and 1.55 fm/c (dotted line), respectively. 

In Fig. 1(a) is shown the pt spectrum of direct pho- 
tons. We can see, if we count thermal photons from 
To = 0.35 fm/c, then the curve is higher than data. 
The rest three curves with delayed counting can meet 
the data points. This is to say, QGP is formed later than 
the hydro initial time and the photon emission rate dur- 
ing [0.35, 0.75] fm/c is overestimated. The measured pt 
spectrum rules out 0.35 fm/c, but can not tell which of 
the rest three is the best beginning time. More precise 
measurement like elliptic flow is needed. 

In Fig. 1(b), prompt photons (black dashed line) and 
thermal photons (curve type same as in Fig. 1(a)) are 
shown separately. All the four thermal curves meet each 
other at low pt region and also meet ALICE data points 
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there. Why? Because the hydro evolution at late stage 
is common for the four curves, which is well controlled 
with bulk hadronsfioj. 

In Fig. 1(d) is shown the elliptic flow of thermal pho- 
tons, which increase from blow to up, for the four begin- 
ning time of counting thermal photons, 0.35, 0.75, 1.15 
and 1.55 fm/c, respectively. This is a clear illustration 
how the elliptic flow grows with time. The initial flow 
velocity is zero, so that the photons emitted at the first 
time interval have zero elliptic flow. Thus the red line 
decreases to zero at high pt region. 

The elliptic flow of direct photons is shown in Fig. 1(c). 
Because of the dominance of prompt photons, it vanishes 
at high Pt , for all the four choices of beginning time. Also 
for the four curves, the increase and decrease mode and 
the turning point at pt = 2 GeV, agrees well with ALICE 
data (full squares). Among the curves, the elliptic flow 
can increase a few times due to the delayed counting of 
thermal photon emission, a much stronger effect than 
fluctuating the initial condition. This big increase also 
tells elliptic flow is more sensitive than pt spectrum to 
the QGP formation time. 

One may expect to meet the large elliptic flow observed 
by ALICE via even later beginning time. This can en- 
large the elliptic flow of thermal photons, but weaken 
the percentage of thermal photons in eq. (4). Thus the 
elliptic flow of direct photons can not be increased. Be- 
sides, too late beginning may also be ruled out by the pt 
spectrum. 

The resulted elliptic flow of direct photons, even with 
delayed 3-flavour QGP photon counting, still below the 
data, tells us the limitation of the smoothed 3+lD ideal 
hydrodynamical model. 

Now we turn to RHIC. In Fig. 2(a) is shown the pt 
spectra of direct photons at centralities 0-20% and 20- 
40%, where the results of 0-20% centrality have been 
increased with a factor of ten for better illusion. The 
PHENIX results are plotted as full triangles and full 
squares for the two centralities, respectively. The solid 
lines correspond to the conventional choice, thermal pho- 
tons begin to emit at tq = 0.6 fm/c. The simplified cal- 
culation, counting prompt photons (without energy loss) 
and thermal photons, does reproduce the PHENIX ob- 
served Pt spectra, as well as the consideration of four 
main contributions with detailed consideration of energy 
loss in [isj . Surprisingly, all the four kinds of curves, for 
four choices of beginning time such as 0.6, 1.1, 1.6 and 
2.1 fm/c, are so close to each other that the measured pt 
can not distinguish them. 

In Fig. 2(b) and (c). is shown the elliptic flow of 
direct photons for centrality 0-20% and 20-40%, re- 
spectively. Similar to the case of PbPb collisions at 
y/sNN =2.76 TeV, later counting of thermal photon 
emission can increase the elliptic flow. Comparing with 
To = 0.6 fm/c, the beginning time of 2.1 fm/c makes the 
enhance of elliptic flow by a factor of about 3 for central- 
ity 0-20% and 4 for centrality 20-40%! This is a much 
bigger effects than the initial fluctuation. 



The pt-dependence of calculated V2 has a similar 
tend as the observed data, increase first, then decrease, 
due to the competition between thermal photons and 
prompt photons. However, the turning point appears 
at pt=l-5 GeV/c, smaller than both PHENIX data and 
the value calculated for PbPb collisions at 2.76 TeV. 

There are two possible reasons for the smaller turn- 
ing Pt- One is that the initial time of AuAu system is 
To=0.6 fm/c, later than tq = 0.35 fm/c in PbPb sys- 
tem. The big tq cuts the earliest growth of the flow ve- 
locity and weakens the competition of thermal photons, 
both in the magnitude of elliptic flow and in the slope 
of Pt spectrum. An early initial time of hydrodynami- 
cal expansion, can certainly improve the situation. The 
other reason may be the equations of state(EoS), the re- 
lation between pressure and energy density, is based on 
first-order transition in AuAu system, while a more re- 
alistic EoS with cross-over phase transition from lattice 
QCD [lOl, in PbPb system. A relatively long co-existence 
of the two phases appears in AuAu system, while no pres- 
sure to accelerate the expansion in the coexistence region. 

Again, though the later thermal photon emission can 
increase the elliptic flow of direct photons, the best choice 
of beginning time, ie. 2.1 fm/c, makes the elliptic flow 
still below the PHENIX data points, consistent with the 
case of PbPb collisions at 2.76TeV. 



IV. CONCLUSION AND DISCUSSION 

In conclusion, we calculated pt spectrum and elliptic 
flow of direct photons from AuAu collisions at 200 GeV 
and PbPb collisions at 2.76 TeV, with a 3+lD event- 
averaged smoothing hydrodynamical description with a 
constrain of big hadronic database, especially the elliptic 
flow of hadrons. The obtained results have been com- 
pared with the recent direct photon data from RHIC and 
LHC. 

The measured large elliptic flow of direct photons (con- 
sistent between RHIC and LHC) is still under-predicted 
with present hydrodynamical model. More details such 
as fluctuation, viscosity and earlier beginning of collective 
motion can improve the prediction. The earlier beginning 
of collective motion, i.e., tq < 0.6 fm/c can improve the 
competition of thermal photons to prompt ones, as we 
can see from the two different initial times in the two hy- 
drodynamical models used for AuAu collisions at RHIC 
energy and PbPb collisions at LHC energy. 

The most important point we learned here is, a delayed 
thermal photon emission from 3-flavour QGP, can explain 
the elliptic flow much better and be allowed by measured 
Pt spectra of direct photons. The best beginning time to 
count photons from 3-flavour QGP is about 2 fm/c, later 
than the beginning of collective motion. 

In our calculation, the pre-equilibrium contribution is 
not calculated separately. The elliptic flow of this con- 
tribution can not be big to solve the large elliptic flow 
puzzle, because of the very early emission time. In our 
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Figure 2: (Color Online) Direct photons from Au+Au collisions at ^snm = 200 GeV. PHENIX data are shown as full triangle 
and full squares for centrality 0-20% and 20-40%, respectively, (a): Transverse momentum spectrum of direct photons, Results 
of centrality 0-20% have been increased with a factor of ten for better illustration. PHENIX data from (b) and (c): 
Compare calculated elliptic flow 112 of direct photons with measured dataQ for centrality 0-20% and 20-40%, respectively. 



calculation, the beginning time of counting thermal pho- 
tons is adjusted to mimic the experimental data, which 
cover the pre-equilibrium contribution. But the emis- 
sion rate is overestimated when the equilibrium one is 
taken as an approximation. That's why a break of the 
emission between prompt photons and thermal photons 
appears in the calculation. As a result, this brings an 
error for our estimation of QGP formation time, which 
should be slightly later than the best beginning time 
we have extracted. Deeper study of the pre-equilibrium 
stage, with more knowledge of the matter before equi- 
librium, will improve our estimation of QGP formation 
time. And what happens before the QGP formation in 
the collision system, becomes an even more interesting 
topic. Cross checking, i.e. from the nuclear modifica- 



tion factor and elliptic flow of hadrons at high pt and 
heavy-flavour hadrons, should be done to test this ear- 
lier collective expansion and later formation of 3-falvour 
QGP in high energy nuclear collisions. 
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